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MILITARY  CONSTRUCTION  ENGINEERING 
AND  DESIGN  COST  FORECASTS 


1 INTRODUCTION 
Background 

Each  year,  the  Directorate  of  Military  Construc- 
tion (DMC)  of  the  Office  of  the  Chief  of  Engineers 
(OCE)  establishes  annual  targets  for  each  Corps  of 
Engineers  Division/District’s  engineering  and  design 
(E&D)  costs  for  military  construction.  These  targets, 
which  are  expressed  as  a percent  of  the  dollar  value 
of  construction  designed,  are  currently  established 
by  an  empirical  procedure  based  on  the  previous  4 to 
5 years  of  Division /District  performance  and  the 
estimated  cost  of  construction  (ECC)  for  the  next 
fiscal  year. 

Objective 

The  objective  of  this  work  is  to  develop  a statistical 
model  for  forecasting  E&D  costs  to  aid  DMC  in 
establishing  E&D  targets  for  each  of  its  Divisions 
and  Districts.  The  model  is  to  be  formulated  so  that 
easily  obtainable  data  can  be  used. 

Approach 

Data  available  on  the  variables  thought  to  affect 
E&D  costs  were  reviewed.  A model  was  formulated 
based  on  these  variables,  and  regression  analysis  was 
used  to  fit  the  model  to  the  data.  The  fit  of  the  model 
was  checked,  and  its  capability  to  accurately  forecast 
E&D  rates  was  determined. 


4 METHOD  OF  DATA  REVIEW 

E&D  costs  were  collected  from  ti  e OCE  Program 
Review  and  Analysis:  Division  and  District  Perform- 
ance Data  report  for  the  10-year  period  FY  66 
through  FY  75.  Appendix  A presents  the  raw  data. 
Prior  to  FY  71.  the  estimated  cost  of  construction 
was  defined  as  85  percent  of  the  programmed  cost, 
or  if  the  programmed  cost  was  not  established,  the 
equivalent  cost.  In  the  second  quarter  of  FY  71,  this 
percentage  was  changed  to  90  percent.  The  raw  data 
were  adjusted  to  reflect  this  change,  as  shown  in 
Appendix  A. 

The  following  Divisions/Districts  were  eliminated 
from  consideration  because  10  years  of  data  were  not 


available  for  them: 

a.  European  Division 

b.  Huntsville  Division 

c.  Kansas  City  District 

d.  Norfolk  District 

e.  Los  Angeles  District 

The  Mediterranean  Division  was  eliminated  be- 
cause the  recent  Saudi  Arabian  Government  work- 
load has  invalidated  its  historical  data.  The  follow- 
ing Divisions/ Districts  had  sufficient  data  for  analy- 


a.  Baltimore  District 

b.  Fort  Worth  District 

c.  New  York  District 

d.  Omaha  District 

e.  Alaska  District 

f.  Mobile  District 

g.  Savannah  District 

h.  Sacramento  District 

i.  Pacific  Ocean  Division 


O MODEL  FORMULATION  AND  ANALYSIS 

Based  on  the  initial  data  review  and  the  detailed 
formulation  procedure  shown  in  Appendix  B,  the 
relationship  between  E&D  costs  for  a particular 
Division/District,  estimated  cost  of  construction, 
and  time  was  postulated  as: 

Dj  - bf}i  + bIiC  + b2iT+b3iTC  [Eq  1] 

where  Dj  = predicted  E&D  costs  for  the  rth  Divi- 
sion/District ($  MIL) 

C = estimated  cost  of  construction  ($  MIL) 
T = time  period  (FY  66  =0,  FY  67  = 1 . . .) 

7C  = 7xC 

bQj,  bjj.  b2i,  b3i  = coefficients  for  the  rth 
Division/District. 

The  E&D  rate  (P-)  is  defined  as  100  x D/C. 

Because  of  the  limited  amount  of  data  available 
for  individual  Divisions/Districts,  a test  was  made  to 
determine  whether  the  data  from  some  or  all  of  the 
Districts  could  be  combined  because  of  similar  E&D 
cost  performance.  (Appendix  B details  the  test  pro- 
cedure.) It  was  found  that  the  data  from  five  Dis- 
tricts (Alaska,  Baltimore,  Fort  Worth,  New  York, 
and  Omaha)  could  be  pooled,  thereby  increasing 
confidence  in  the  prediction  equations.  The  regres- 
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sion  analyses  for  the  remaining  Divisions/Districts 
were  based  only  on  their  own  data. 

The  postulated  model,  when  regressed  against 
FY  66  through  FY  77  data,  provided  a good  fit  for 
eight  Districts,  explaining  92.4  to  99.8  percent  of  the 
variance  of  the  original  data  (Table  1).  The  model 
did  not  fit  the  data  for  Pacific  Ocean  Division;  it 
explained  only  44.8  percent  of  the  variance  of  the 
data  and  hence  is  considered  inadequate.  The  large 


variability  and  erratic  behavior  in  that  Division’s 
E&D  cost  data  caused  the  lack  of  fit. 

4 MODEL  RESULTS  AND  VERIFICATION 
Results 

Figures  1 through  8 show  the  FY  76  and  FY  77 
prediction  equations  and  graphs  for  the  eight  Dis- 


Tabtc  1 

Ryiwlod  Retail* 

D=/>,  + 6,C  + 6,r 


Dh/Dbtrict 

FYData 

bx 

bx 

Standard  Error 
of  Eatln  ate, 
S(SMIL) 

Percent  of  Variance 
of  Original  Data 
Explained  by  Model 

Rl 

Alaska 

1966-1972 

-.0215 

.0490 

.0688 

.160 

.997 

1966-1973 

.3605 

.0410 

.0048 

.153 

.998 

1966-1974 

.2963 

.0424 

.0143 

.150 

.998 

1966-1975 

.3483 

.0413 

.0062 

.149 

.998 

Baltimore 

1966-1972 

.1826 

.0465 

.0908 

.160 

.997 

1966-1973 

.1399 

.0474 

.0906 

.163 

.998 

1966-1974 

.1601 

.0479 

.0690 

.150 

.998 

1966-1975 

.1835 

.0489 

.0330 

.149 

.998 

Fort  Worth 

1966-1972 

.6371 

.0377 

.2537 

.160 

.997 

1966-1973 

.7278 

.0363 

.2540 

.153 

.998 

1966-1974 

.7178 

.0364 

.2537 

.150 

.998 

1966-1975 

.8602 

.0337 

.2733 

.149 

.998 

Mobile 

1966-1972 

.3067 

.0413 

.248 

.935 

1966-1973 

.1165 

.0466 

.599 

.732 

1966-1974 

-.1259 

.0504 

.564 

.870 

1966-1975 

-.2325 

.0520 

.532 

.924 

New  York 

1966-1972 

.6895 

.0369 

.0545 

.160 

.997 

1966-1973 

.5357 

.0398 

.0679 

.153 

.998 

1966-1974 

.5545 

.0397 

.0621 

.150 

.998 

1966-1975 

.5772 

.0397 

.0531 

.149 

.998 

Omaha 

1966-1972 

.1257 

.0462 

.1516 

.160 

.997 

1966-1973 

.1327 

.0460 

.1550 

.153 

.998 

1966-1974 

.0827 

.0465 

.1631 

.150 

.998 

1966-1975 

.0675 

.0461 

.1779 

.149 

.993 

Pacific  Ocean  Division 

1966-1972 

2.9405 

.0245 

1.686 

.271 

1966-1973 

2.2439 

.0300 

1.584 

.433 

1966-1974 

1.7340 

.0338 

1.540 

.499 

1966-1975 

1.9040 

.0306 

1.510 

.448 

Sacramento 

1966-1972 

-.0693 

.0622 

.228 

.969 

1966-1973 

.0581 

.0588 

.235 

.972 

1966-1974 

.1485 

.0565 

.235 

.975 

1966-1975 

.2119 

.0551 

.226 

.985 

Savannah 

1966-1972 

1.3901 

.0244 

.368 

.359 

1966-1973 

1.3361 

.0256 

.342 

.383 

1966-1974 

1.3962 

.0246 

.317 

.662 

1966-1975 

.0100 

.0459 

.497 

.925 
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tricts  for  which  the  model  provided  a good  fit.  The 
prediction  limits  indicate  the  accuracy  of  the  predic- 
tion. The  following  example  illustrates  how  they  are 
interpreted.  For  Baltimore  District  (Figure  2),  if  the 
estimator  assumes  an  estimated  cost  of  construction 
of  S160  million  in  FY  76,  he/she  can  be  95  percent 
confident  that  the  actual  E&D  cost  will  be  between 
$8.0  and  $8.8  million.  Equivalently,  he/she  can  be 
95  percent  confident  that  the  actual  E&D  rate  will  be 
between  5.0  and  5.5  percent. 

The  graphs  have  been  plotted  over  a reasonable 
range  of  estimated  construction  costs  for  each  Dis- 
trict, since  excessive  extrapolation  is  risky.  For 
example,  New  York  District,  which  has  been  design- 
ing an  average  of  $40  to  $45  million  of  construction 
per  year,  would  not  be  expected  to  perform  in  the 
same  manner  if  it  were  suddenly  tasked  to  design 
$200  million  of  construction  per  year. 

Verification 

The  validity  of  a forecasting  model  is  determined 
by  its  capability  to  predict  future  performance.  A 
systematic  approach  to  verifying  the  model  is  to  con- 
duct a retrospective  test.  Since  lack  of  the  current 
year’s  data  when  the  next  year's  E&D  limits  must  be 
established  may  require  using  the  ryiodel  to  predict 


performance  2 years  in  the  future,  the  model  was 
tested  for  predicting  both  1 and  2 years  ahead.  The 
retrospective  test  was  conducted  by  comparing  the 
actual  FY  74  E&D  costs  to  the  FY  74  costs  predicted 
by  the  model  developed  from  FY  66  through  FY  72 
data  (a  2-year-ab*ad  prediction),  and  to  the  FY  74 
costs  predicted  by  the  model  developed  from  FY  66 
through  FY  73  data  (a  1-year-ahead  prediction). 
Analogous  comparisons  were  also  made  between 
actual  FY  75  E&D  costs  and  the  1-  and  2-year-ahead 
predictions. 

Table  2 presents  the  results  of  the  retrospective 
test.  In  all  cases  but  one,  the  actual  E&D  cost  was 
well  within  the  95  percent  confidence  limits.  Savan- 
nah District’s  actual  E&D  cost  for  FY  75  was  above 
the  upper  prediction  limit.  This  result  is  within 
expectations  because  95  percent  confidence  means 
that  approximately  95  out  of  100  actual  E&D  costs 
will  fall  within  the  upper  and  lower  prediction  limits. 

5 MODEL  USE  AND  MAINTENANCE 
Model  Use 

The  model  can  be  used  to  predict  the  performance 
of  a particular  Division/District  based  on  its  past 


Table  2 

Retroepectlve  Teet  Reiult* 


1-Year-Ahead  Prediction  2-Year-Ahead  Prediction 


Actual 

Predicted 

Lower 

Upper 

Predicted 

Lower 

Upper 

E&D 

E&D 

Prediction 

Prediction 

E&D 

Prediction 

Prediction 

Dlv/Dlslrlct 

FY 

Co*t 

Coit 

Limit 

Limit 

Coat 

Level 

Level 

Alaska 

1974 

I 681 

1.601 

1.064 

2.138 

1.965 

1.035 

2.8% 

1975 

1.444 

1.551 

1.139 

1.963 

1.493 

.917 

2.070 

Baltimore 

1974 

7 932 

8.089 

7.648 

8,531 

8.008 

7.516 

8.500 

1975 

6.594 

6.906 

6.460 

7.352 

7.017 

6.462 

7.572 

Fort  Worth 

1974 

8.055 

8.037 

7.563 

8,511 

8.145 

7.460 

8.830 

1975 

10.141 

10.467 

9.882 

11.051 

10.448 

9.663 

11.232 

Mobile 

1974 

6.591 

6.174 

3.881 

8.467 

5.668 

4.642 

6.695 

1975 

7.495 

7.236 

5.384 

9.087 

6.928 

4.283 

9.574 

New  York 

1974 

2.735 

2.822 

2.422 

3.222 

2.739 

2.164 

3.315 

1975 

2.618 

2.783 

2.405 

3.161 

2.822 

2.394 

3 250 

Omaha 

1974 

12.386 

12.050 

11.645 

12.455 

12.077 

11.529 

12.625 

1975 

11  422 

11.314 

10.910 

11.719 

11.185 

10.740 

1 1 .630 

Sacramento 

1974 

: 06S 

5.370 

4.633 

6.107 

5.558 

4.703 

6.414 

1975 

6 72 i 

6 924 

6.132 

7.716 

7.110 

6 174 

8.046 

S.i'annah 

1974 

3.970 

4.029 

2.392 

5.66t 

3.959 

2 052 

5.866 

1975 

8 266 

5 567 

3.779 

7 355 

5.674 

2.101 

9.248 
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performance.  DMC  must  exercise  its  managerial 
control  in  establishing  E&D  targets  which  may  or 
may  not  be  equal  to  the  predicted  rate.  Since  District 
performance  is  influenced  by  the  established  target, 
DMC  may  predict  the  expected  E&D  rate  and  then 
establish  the  target  either  below  or  above  this  value 
to  restrict  or  expand  the  E&D  effort,  as  desired.  For 
example,  if  restrictions  are  deemed  necessary  for  a 
particular  District  based  on  past  performance,  the 
DMC  may  wish  to  establish  the  target  at  the  lower  95 
percent  confidence  limit. 


parameter  changes  as  additional  data  are  added, 
beginning  with  the  data  through  FY  72.  The  coeffi- 
cient for  the  estimated  cost  of  construction  remains 
fairly  stable;  however,  in  some  Districts,  such  as 
Alaska,  the  other  parameters  change  as  additional 
data  are  added.  This  indicates  that  although  the 
basic  model  may  not  be  changing,  the  effects  of  the 
variables — particularly  time — on  E&D  costs  may  be. 
and  caution  should  be  exercised  in  predicting  several 
years  in  advance.  New  parameters  should  be  deter- 
mined when  each  year’s  data  are  added. 


The  prediction  equations  given  for  each  of  the 
Districts  in  Table  1 (FY  66  through  FY  75  data)  can 
be  used  to  calculate  an  estimated  E&D  cost  by  sub- 
stituting the  estimated  cost  of  construction  and  the 
appropriate  time  period  (FY  76  = 10,  FY  77  » 11, 
etc.).  The  E&D  rate  is  calculated  by  dividing  the 
E&D  cost  by  the  estimated  cost  of  construction.  In 
Figures  1 through  8,  the  value  of  time  is  already 
accounted  for  and  combined  with  the  constant  term; 
therefore,  only  an  estimated  cost  of  construction  is 
required.  For  example,  Baltimore’s  prediction  equa- 
tion is: 


D = .1835  + .0489  (O  + .0330  (7) 

Assuming  estimated  construction  costs  for  FY  76  at 
$160  million, 

D = .1835  + .0489  (160)*  + .0330  (10).** 


The  goodness  of  fit  of  the  model  can  be  checked  by 
examining  the  estimated  standard  error  of  the  esti- 
mate, S,  and  the  percent  of  variance  in  the  original 
data  explained  by  the  model,  R1  (Table  1).  An  in- 
creasing value  of  S or  a decreasing  value  of  R 2 would 
indicate  that  the  basic  model  itseli  may  be  changing 
and  steps  should  be  taken  to  determine  where 
changes  are  occurring.  Stable  or  decreasing  trends  in 
S indicate  a stable  model.  The  high  values  for  R1 
indicate  the  model  is  describing  the  data  very  well. 

Finally,  the  accuracy  of  the  model  as  a predictor 
should  be  checked  annually.  Actual  values  which 
continually  fall  outside  the  confidence  limits  may 
indicate  changes  in  the  model  itself. 
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CONCLUSIONS 

AND  RECOMMENDATIONS 


Therefore, 


D = $8,337  million 


and 


The  model  presented  in  this  report  is  a valid  tech- 
nique for  forecasting  future  E&D  costs.  Predictions 
were  presented  for  FY  76/FY  77  for  the  eight  Mili- 
tary Construction  Divisions /Districts  which  had 
sufficient  data  for  analysis. 


P = ~^-xi°°=5.2  percent. 

Based  on  this  value  of  expected  actual  performance, 
the  established  limit  may  be  set  according  to  the  cur- 
rent management  policy  for  the  District. 

Model  Maintenance 

The  stability  of  the  model  parameters,  the  good  fit 
of  the  model  with  the  data,  and  the  accuracy  with 
which  the  model  is  predicting  actual  E&D  costs 
should  be  checked  periodically.  Table  1 shows  the 


*160  is  used  because  the  model  ts  based  on  millions  of  dollars. 

**FY  7(,  = |0. 


It  is  recommended  that  the  DMC  use  this  model 
as  an  aid  in  establishing  E&D  targets.  Although  the 
model  can  be  used  to  predict  costs  1 or  2 years  in 
advance,  the  best  results  are  obtained  for  predictions 
made  1 year  ahead.  Thus,  the  model  should  be  up- 
dated with  the  past  year’s  data  before  determining 
the  new  year’s  targets. 

Because  the  E&D  costs  predicted  by  the  model  are 
based  on  past  Division/District  performance  and 
conditions,  significant  changes  in  performance  or 
conditions  could  invalidate  the  historical  data  on 
which  the  model  is  based.  Therefore,  the  DMC  must 
consider  changes  in  operating  conditions  or  manage- 
ment policy  within  the  Division/District  when  estab- 
lishing E&D  targets. 
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Figure  la.  Prediction  of  E&D  costs/rates — Alaska — FY  76. 


Figure  lb.  Prediction  of  E&D  costs/iates — Alaska — FY  77 


ws&rvrzr.  ^T^7,  - -» 


(90 1 * $ ) 1S00  0 Puo  3 


E end  D COST  vs  E and  D RATE  vs 

ESTIMATED  COST  OF  CONSTRUCTION  ESTIMATED  COST  OF  CONSTRUCTION 


Figure  3b.  Prediction  of  E&D  costs/rates — Fort  Worth — FY  77. 


E and  D COST  vs  E and  D RATE  vs 

ESTIMATED  COST  OF  CONSTRUCTION  ESTIMATED  COST  OF  CONSTRUCTION 


Prediction  of  E&D  costs/rates — New  York — FY  76. 


E and  D COST  vs  E and  D RATE  vs 


Figure  6b.  Prediction  of  E&D  costs/rates — Omaha — FY  77. 


E and  D COST  vs  E and  D RATE  vs 


Figure  8.  Prediction  of  E&D  costs/rates — Savannah — FY  76/FY  77. 
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FORT  WORTH  DISTRICT 


FY 

ECC 
($  Mlt) 

AdJ  ECC 
(SMU) 

E&DCoat 

(SMU) 

EADRale 

<%) 

1966 

83.034 

66.742 

3.285 

4.922 

1967 

56.627 

59.958 

2.960 

4.937 

1988 

52.366 

55.446 

3.16C 

5.714 

1%9 

74.778 

79.176 

4.435 

5.601 

1970 

53.133 

56.258 

3.969 

7.055 

1971 

75.350 

75.350 

4.598 

6.102 

1972 

104.900 

104.900 

6.128 

S.842 

1973 

122.985 

122.985 

6.921 

5.628 

1974 

145.416 

145.426 

8.055 

5.539 

1975 

204.842 

204.842 

10.141 

4.951 

HUNTSVILLE  DIVISION \ 


FY 

ECC 

(SMtl) 

A«U  ECC 
(SMU) 

EttDCoat 

(SMU) 

EiD  Rate 

(%) 

1966 

1967 

1968 

1969 

1970 

204.532 

216.563 

8.184 

3.779 

1971 

241.000 

241.000 

9.773 

4.055 

1972 

216.700 

216.700 

9.369 

4.323 

1973 

91.836 

91.836 

4.875 

5.308 

1974 

48.527 

48.527 

2.426 

4.999 

1975 

180.712 

180.712 

11.041 

6.110 

KANSAS  CITY  DISTRICT 


FY 

ECC 
(S  MU) 

AdJ  ECC 
(S  MID 

E&DCo*t 

(SMU) 

E&DRate 

<%) 

1966 

22.335 

23.649 

931 

3.937 

1967 

27.840 

29.478 

1.300 

4.410 

1968 

62.094 

65.746 

1.939 

2.949 

1969 

50.116 

53.064 

2.169 

4.088 

1970 

1971 

1972 
197.1 
1974 

20.920 

22.151 

1.276 

5.761 

5 391 


5.843 


5.84.1 


.115 


fjc 


v 

£ 

£ 
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LOS  ANGELES  DISTRICT 


FY 

ECC 
(S  Mil) 

AdJ  ECC 
(SMU) 

E&D  Coat 
(SMU) 

EtDRdt 

(%) 

1966 

33.938 

35.934 

1.910 

5.315 

l%7 

40.013 

42.367 

2.284 

5.391 

l%8 

44.800 

47.435 

2.65? 

5.591 

1%*) 

57.364 

60.738 

4.039 

6.650 

l<T0 

33.373 

35.336 

2.336 

6.611 

1971 

3.257 

3.257 

.318 

9.764 

1972 

1973 

1974 

1975 

.658 

.658 

.110 

16.717 

NEW  YORK  DISTRICT 


FY 

ECC 

(SMU) 

AdtJECC 

(SMB) 

EADCoat 

(SMI) 

EAD  Rate 

(%) 

1966 

52.190 

55.260 

2.723 

4.928 

1967 

47.350 

50  135 

2.549 

5.084 

1968 

34.493 

.36.522 

2.035 

5.572 

1969 

33.334 

35.295 

2.259 

6.400 

1970 

30.393 

32.181 

2.382 

7.402 

1971 

27.104 

27.104 

1.735 

6.401 

1972 

44.376 

44.376 

2.638 

5.945 

1973 

60.947 

60.947 

3.470 

5.693 

1974 

43.805 

43.805 

2.735 

6.244 

1975 

42.084 

42.084 

2.618 

6.221 

{ 

! 


MEDITERRANEAN  DIVISION 

NORFOLK  DISTRICT 

FY 

ECC 

($  Mil) 

AdJ  ECC 
(S  MU) 

E&D  Co«t 
(SMU) 

E&D  Rate 

(%) 

FY 

ECC 

(SMU) 

Adj  ECC 

(SMU) 

E&D  Coat 
(SMU) 

E&D  Rate 

(%) 

1966 

42.518 

45.019 

1.727 

3.836 

1966 

31.611 

33.470 

1.739 

5.1% 

19(>7 

81  305 

86.087 

3.490 

4.054 

1967 

42.120 

44.597 

2.361 

5.294 

1968 

49.106 

51.994 

2.076 

3.993 

1968 

53.869 

57.038 

2.968 

5.204 

1969 

36.544 

35.694 

1.248 

3.225 

1969 

28.139 

29.794 

1.993 

6.689 

1970 

20.820 

22.045 

.907 

4.114 

1970 

29.068 

30.778 

2.306 

7.492 

1971 

13.844 

18.844 

.988 

7.187 

1971 

1972 

27.445 

27.445 

1.871 

6.817 

1972 

1973 

20.695 

20.695 

1.646 

7.954 

1973 

1974 

96.356 

96.356 

3.165 

3.285 

1974 

69.945 

69.945 

2.730 

3.903 

1975 

1259  372 

1259  372 

9.505 

.755 

1975 

76.553 

76.553 

3.522 

4.601 

MOBILE  DISTRICT 

OMAHA  DISTRICT 

FY 

ECC 
(S  MU) 

Adj  ECC 
(SMU) 

E&D  Coil 
(SMU) 

E&D  Rate 

(%) 

FY 

ECC 

(S  Mil) 

Adj  ECC 
($  MU) 

E&D  Coat 
(SMU) 

E&D  Rate 

<%) 

l%(> 

72  269 

76  520 

3.206 

4.190 

1966 

37.368 

39.778 

2.085 

5.242 

1%' 

'9  623 

04.518 

4.219 

4.992 

l%7 

29.604 

31.345 

1.572 

5.015 

1968 

49  803 

52.732 

2.370 

4.494 

l%8 

44.580 

47.202 

2.560 

5 423 

1969 

30  889 

32  706 

1 821 

5 568 

l%9 

28  811 

29.976 

2.075 

6 922 

|9'0 

63.821 

67.575 

2.976 

4,404 

1970 

72.893 

77.181 

4.211 

5.456 

I9-| 

7 1 . 26 1 

71  261 

3.162 

4.437 

I97| 

123.120 

123.120 

6.697 

5.439 

19' 2 

95  799 

95.799 

4.241 

4 427 

1972 

167  840 

167  840 

8.747 

5.212 

I9'l 

79.929 

'9.929 

5.087 

6.364 

I973 

228  1 ’5 

228  175 

11.694 

5.125 

|9'4 

129  %4 

129  %4 

6.591 

5.071 

I974 

232.203 

232.203 

12.386 

5.334 

19-5 

146.148 

146  148 

7.495 

5.128 

19"5 

210.027 

210.027 

11.422 

5.438 
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PACIFIC  OCEAN  DIVISION 


SAVANNAH  DISTRICT 


FY 

ECC 

<SMB 

A4JECC 

(SMB 

EADCorf 

(SMB 

EAD  Rate 

(%) 

FY 

ECC 

(SMB 

Ac#  ECC 
(SMB 

EADCwt 

(SMB 

EAD  Rate 

<%) 

l%6 

143.306 

151.735 

5.275 

3.476 

1966 

67.693 

71.675 

2.723 

3.799 

l%7 

I39.«  55 

147.372 

5.423 

3.680 

1967 

67.615 

71.592 

2.886 

4.031 

1968 

66.728 

93.947 

4.944 

5.263 

1968 

55.518 

58.784 

2.745 

4.670 

1969 

135.920 

143.915 

9.132 

6.345 

1969 

58.247 

61.673 

3.005 

4.872 

1970 

73.209 

77.515 

6.490 

8.373 

1970 

59.922 

63.447 

3.489 

5.499 

1971 

(>6.213 

66.213 

3.871 

5.846 

1971 

41.497 

41.497 

2.201 

5.304 

1972 

74.332 

74.332 

3.970 

5.341 

1972 

66.512 

66.512 

3.2  5 

4.954 

197.1 

10.445 

40.44S 

2.743 

6.782 

1973 

67.684 

67.684 

3.204 

4.734 

1974 

56.003 

56.003 

2.S27 

4.512 

1974 

105.320 

105.320 

3.970 

3.769 

1975 

133.867 

133.867 

4.844 

3.619 

1975 

169.654 

169.654 

8.266 

4.872 

SACRAMENTO  DISTRICT 


FY 

ECC 

(SM11) 

A«y  ECC 
(SMB 

EADCoct 

(SMB 

EADRa 

(%) 

1966 

25.684 

27.195 

1.609 

5.917 

1967 

23.844 

25.247 

5.652 

1968 

46.387 

49.115 

2.637 

5.369 

1969 

33.090 

35.036 

2.248 

6.416 

1970 

34.985 

37.043 

2.453 

6.622 

1971 

71.446 

71.446 

4.598 

6.436 

1972 

68.049 

68.049 

4.035 

5.930 

1973 

81.456 

81.456 

4.643 

5.700 

1974 

90.403 

90.403 

5.068 

5.606 

1975 

120.013 

120.013 

6.721 

5.600 

wwffft*  vv«  ,<«n 


APPENDIX  B: 


Figure  B2  is  a graphical  interpretation  of  Eq  B3. 


MATHEMATICAL  DEVELOPMENT 
Model  Formulation 

Two  variables  were  available  for  use  in  estimating 
E&D  costs.  D,  for  a particular  Division/District: 
(1)  estimated  cost  of  construction,  C,  and  (2)  time, 
T.  The  first  step  was  to  determine  the  relationship 
between  D,  C,  and  T.  A correlation  matrix  (Table 
Bl)  was  computed  from  the  adjusted  data. 

Table  Bl 

Correlation  of  Depeodait/In  depen  dent  Variable* 


D 

C 

r 

D 

1.000 

.969 

.482 

C 

.969 

1.000 

.437 

T 

.482 

.437 

1.000 

The  initial  Division /District  model  at  a base  time 
was  formulated  as  a linear  relationship  between  D 
and  C: 


Figure  B2.  E&D  rate  vs  estimated  cost 
of  construction. 

The  average  percent  E&D  decreases  as  the  esti- 
mated cost  of  construction  increases.  This  phenome- 
non, which  is  known  as  ‘'economy  of  scale,”  is  due  to 
the  amortization  of  the  fixed  charge  over  a larger 
workload  base.  If  there  is  no  economy  of  scale  within 
a Division/District,  the  general  model  would  auto- 
matically collapse  to  Eq  B4  and  BS  and  Figure  B3. 


D =AjC 

F - 100  a, 


[Eq  B4] 
[Eq  B5] 


D-aa  + a,C 


(Eq  Bl] 


The  E&D  costs  equal  some  fixed  costs,  aQ,  plus  some 
fixed  percentage,  av  of  the  estimated  cost  of  con- 
struction (Figure  Bl). 


Figure  Bl.  E&D  costs  vs  estimated  cost 
of  construction. 

The  average  E&D  rate  (P  ) is  defined  as  the  E&D 
cost  divided  by  the  estimated  cost  of  construction 
designed  times  100: 

IEqB2| 

The  model  can  then  be  interpreted  in  terms  of  aver- 
age percentages  bv  dividing  Eq  Bl  by  C: 

F = 100  al  + WO  df/C  [Eq  B3] 


The  next  variable  to  bt  r ansidered  was  time.  Costs 
of  both  design  and  construction  increase  over  time, 
probably  at  different  rates.  Some  of  these  costs  are: 

a.  Materials 

b.  Labor 

c.  Effects  of  Environmental  Impact  Assessments/ 
Statements 

d.  Workload  levels 

e.  Manpower  levels. 

To  account  for  changes  in  D and  C across  time  the 
total  change  in  each  was  modeled  by  uniform  annual 
rates  compounded  over  time.  The  advantage  of  this 
method  is  that  it  does  not  requite  the  imposition  of 
indices  that  may  not  completely  correspond  to  the 
inflationary  pressures  exerted  on  the  type  of  design 
and  construction  that  the  Corps  performs.  These 
uniform  annual  rates  reflect  not  only  inflationary 
effects  on  D and  C,  but  also  noninflationary  effects 
such  as  changes  in  workload  or  manpower  levels. 

The  uniform  annual  rates  of  change  in  design  and 
construction  costs  can  be  represented  as  ij  and  ic, 
respectively,  where 

Dt  = (/  +i(j){  D t = 0.l .;/«  — / (Eq  B6] 

C,  = (/+ic)'C  t = C 1 „_/  (Eq  B7] 
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Figure  B3.  E&D  costs/rates  vs  estimated  cost  of  construction  without  economy  of  scale. 


The  basic  model  was  generalized  to  reflect  changes 
over  time  as  follows: 

Substituting  Eq  81  into  Eq  B6 

Dt  = (/  + id){  (aQ  + alO 

Solving  Eq  B7  for  C and  substituting 

1 


where  b0  - OqOq 
~ Po°i 


b2  « a^Q 
by  = /?,a, 


Dt  = (/  +.'/{o0  + a,[ 


U+ic)( 

Df-U  + ij)f  u0  + [f+-j-c  1 a‘  Ct 


■]Ct) 


[Eq  B8] 


It  can  be  seen  that  the  general  model  (Eq  B8)  is 
equivalent  to  the  basic  model  (Eq  Bl)  in  the  base 
year  (t  = 0). 

Eq  B8  was  transformed  into  a linear  equation 
using  the  approximation 

(/  ± x)r  —■  uo  ± a\ f 

This  approximation  is  quite  good  for  small  values  of 
X and  when  the  maximum  t is  not  large.  The  follow- 
ing nonlinear  to  linear  transformations  were  made: 

(/  ± id)!  a0  a:  a,t 

1 +id  1 

HTe 

Substituting  the  above  linear  approximations  into 
Eq  E8  gives 

Dt  = {%  ± aj  t)a0  + (p0  * /3(  t)aj  C, 

Dj  = a0ao  + P0OiCf  + flj o0t  T P1a1tCi 
D =b0+b1C  + b1T  + biTC  [Eq  B9] 


and  D and  C are  understood  to  be  costs  in  the  Tth 
year. 

Eq  B9  was  then  generalized  for  the  rth  Division/ 
District  to  the  complete  model  for  E&D  costs. 

Dj  = boi  + bIiC  + b2i T + b3i TC  i = 0,  1. ...  8. 

[Eq  BIG] 

Model  Analysis 

Some  basic  goals  were  established  before  any  com- 
puter regression  runs  were  made:1 

a.  The  final  equation(s)  should  explain  more  than 
90  percent  of  the  variation  (R2  > 0.90) 

b.  The  standard  error  of  the  estimate  should  be 
less  than  10  percent  of  the  mean  (coefficient  of  vari- 
ability < 10  percent) 

c.  All  estimated  coefficients  should  be  statistically 
significant  with  a = 0.05 

d.  There  should  be  no  discernible  pattern  in  the 
residuals. 

The  Statistical  Package  for  the  Social  Sciences 
(SPSS)2  was  used  to  perform  the  regression  analysis. 

The  first  step  in  the  analysis  was  to  determine 
which,  if  any,  of  the  Division/District  models  could 


*N  Draper  and  H.  Smith.  Applied  Regression  Analysis  (John 
Wiley  and  Sons.  Inc..  I%h).  p 238 
2 Norman  H . Nie  cl  al . . SPSS— Statistical  Package  lor  the  Social 
Sciences  (McGraw  Hill.  1975). 
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be  combined  either  in  part  or  wholly  without  affect- 
ing the  statistical  significance  of  the  model.  A modi- 
fied “backward  elimination'*  approach  composed  of 
the  following  steps  was  used:3 

a.  Regress  the  data  against  the  model  with  all  of 
the  variables  included 

b.  Use  the  F-statistic  to  determine  whether  the 
deletion  of  each  variable  would  detract  significantly 
from  the  model 

c.  Use  Bartlett's  Chi-square  statistic  to  determine 
whether  there  is  a significant  difference  between  the 
residual  variances  of  the  districts 

d.  if  residuals  are  homogeneous,  the  procedure  is 
complete.  If  not,  remove  the  data  of  the  Division/ 
District  whose  residuals  are  furthest  from  the  resid- 
ual mean  variance  and  go  to  step  a. 

To  accomplish  the  backward  elimination,  a 
dummy  variable  was  used  to  indicate  from  which 
Division /District  a particular  data  point  came. 
Thus,  the  nine  Divisions/Districts  have  a total  of  36 
possible  variables,  as  shown  in  Eq  BIO. 

The  SPSS  method  required  that  a base  District  be 
chosen  (and  numbered  zero);  the  coefficient  for  the 
other  Districts  represents  their  deviation  from  the 
base  District.  When  the  coefficients  for  all  of  the 
variables  were  determined  through  the  regression 
analysis,  the  coefficients  for  any  given  district  were 
added  to  the  coefficients  for  the  base  District  to  esti- 
mate the  total  coefficient  for  the  given  District. 

Confidence  limits  on  foiecasts  are  a function  of 
the  standard  error  of  the  estimate  and  the  degrees  of 
freedom  related  to  that  estimate.  Degrees  of  freedom 
are  in  turn  based  on  the  number  of  observations  (in 
this  case,  years  of  data).  Each  Division/District  had 

’Draper  and  Smith,  p 167. 


only  10  observations.  The  backward  elimination  pro- 
cedure was  used  to  determine  whether  some  Divi- 
sions/Districts behaved  in  generally  the  same  way 
with  respect  to  E&D  costs.  If  so,  their  data  could  be 
pooled  to  increase  the  degrees  of  freedom,  thus 
narrowing  the  distance  between  the  prediction  limits 
on  any  forecast. 

The  initial  steps  of  the  backward  elimination  pro- 
cedure were  performed  as  follows.  There  were  ini- 
tially 36  possible  regression  coefficients  for  the  nine 
districts  and  their  90  data  points,  as  shown  in  Eq 
BIO.  An  enumeration  tree  (Figure  B4)  was  estab- 
lished to  explicitly  show  the  removal  of  a set  of 
regression  coefficients.  Each  branch  in  the  tree  signi- 
fies the  removal  of  the  associated  coefficient  from  all 
Districts  except  the  base  District,  The  basic  premise 
for  this  is  that  if  the  data  from  the  Districts  are  to  be 
pooled,  all  Districts  are  performing  in  the  same 
manner  and  thus  the  value  of  the  coefficient  for  the 
base  District  (bj0)  adequately  reflects  that  perform- 
ance. 

Table  B2 

Backward  Elimination  Retails 


Model  (I  = 1,8) 

Rl 

df 

F 

a 

1.  All 

.97247 

54 

. 

2.  Without  b0j 

.96192 

62 

2.587 

.02 

3.  Without  bjj 

.96346 

62 

2.209 

.04 

4.  Without 

.96412 

62 

2.047 

.06 

5.  Without  b^j 

.96733 

62 

1.260 

.28 

6.  Without  hj(),  by 

.96589 

63 

1.434 

.20 

The  regression  analysis  was  initially  performed 
with  all  36  regression  coefficients  included  to  provide 
the  best  possible  model  fit,  The  value  of  R1  was  used 
to  determine  goodness  of  fit.  Line  1 of  Table  B2  pre- 
sents the  results  of  this  run.  Next,  runs  of  the  model 
were  made  corresponding  to  each  branch  of  the  tree 
(Figure  B5);  i.e.,  one  set  of  coefficients  was  removed 
at  a time.  Lines  2 through  5 of  Table  B2  show  these 
results.  The  F-statistic  was  calculated  as  shown  in 
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Figure  B4.  Removal  tree  for  one  variable  at  a time. 


Eq  Bi  1 to  determine  whether  there  was  a significant 
decrease  in  the  model  fit  due  to  removal  of  the  set  of 
terms. 

„ _ {Rall  ~ R without^  ^without  ~ ^alP 

' ^l^alMall 


(Eq  Bill 

The  value  of  a represents  the  risk  incurred  in  reject- 
ing the  hypothesis  that  the  removed  terms  are  un- 
necessary, i.e.,  in  concluding  that  the  removed  terms 
were  necessary.  A critical  value  of  a = .10  was  used. 
According  to  the  values  obtained,  only  the  co- 
efficient corresponding  to  the  interaction  term,  can 
be  removed  without  significantly  reducing  the  model 
fit,  (tr  = .28). 

The  next  step  was  to  determine  whether  any  addi- 
tional terms  could  be  removed.  Since  only  the  bjj 
terms  could  be  initially  removed,  bj0  is  the  only  re- 
maining possibility.  Note  that  coefficients  for  the 
base  District  cannot  be  removed  until  the  corre- 
sponding coefficients  for  all  other  districts  have  been 
removed.  Figure  B5  shows  the  removal  of  bjQ\  the 
results  are  shown  in  Table  B2,  line  6.  The  value  of 
a = .20  indicates  no  significant  reduction  in  model 
fit.  For  the  uata  from  the  nine  Districts  under  con- 
sideration, all  of  the  interaction  coefficients  were  the 
same  and  were  assumed  to  be  zero.  The  other  nodes 
in  Figure  B5  have  been  fathomed  (denoted  by  the 
underscoring)  to  indicate  that  since  these  terms 
could  not  be  removed  individually,  they  could 
obviously  not  be  removed  with  other  terms. 


After  the  simplest  model  which  would  not  signifi- 
cantly decrease  the  model  fit  for  the  nine  Districts 
was  determined,  determining  whether  the  data 
could,  in  fact,  be  pooled  was  the  next  step.  This  was 
done  using  Bartlett’s  Chi  square4  test  to  determine 
whether  the  residual  variance  differed  from  District 
to  District.  The  variances  of  the  residuals  for  each 
Division/District  are  shown  in  Table  B3  (Step  1). 
According  to  Bartlett’s  test,  the  District-to-District 
variations  differ  significantly  (probability  = 1.0); 
therefore,  the  data  from  all  nine  Districts  could  not 
be  pooled. 

The  data  for  the  District /Division  with  the  highest 
variation  (in  this  case  Pacific  Ocean)  were  removed 
and  the  process  was  repeated. 

The  results  of  this  process  through  five  iterations 
are  shown  in  Table  B3.  Four  Districts  had  to  be  re- 
moved before  the  data  for  the  remaining  five  could 
be  pooled. 


Table  B3 


Result!  of  Test  for  Homogeneous  Residual  Variance 
Between  Districts/Divisions 


Dlstrlct/Dlvlilon  Step  1 

Step  2 

Step  3 

Step  4 

StepS 

Baltimore 

.01387 

.01387 

.04032 

.01387 

.01387 

Fort  Worth 

.01744 

.01744 

.09565 

.0P44 

.01744 

New  York 

.02078 

.01917 

.023% 

.02078 

.02078 

Omaha 

.01705 

.01705 

.01811 

.01705 

.01705 

Alaska 

.01738 

.01738 

.01907 

.01738 

.01739 

Sacramento 

.0381 1 

.03813 

.03910 

.03812 

Mobile 

.18558 

.18558 

.18568 

Savannah 

.19887 

19887 

Pacific  Ocean 

1.93218 

Probability* 

1 00 

1.00 

1.00 

.30 

.01 

•Probability,  based  on  Bartlett’s  Chi  square  statistic,  that  there 
is  a significant  difference  between  the  variances  from  District  to 
District 


The  same  process  was  used  to  determine  if  any  of 
the  four  removed  Districts  cou'd  be  pooled.  It  was 
found  that  each  of  these  had  to  be  treated  separately. 


4G.  W.  Sncdecor  and  W G Cochran.  Statistical  Methods 
(Iowa  State  University  Press.  1971).  pp  296-298 
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